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Abstract 

Most quantum entanglement investigations are focused on two qubits or some finite (small) chain structure, since the 
infinite chain structure is a considerably cumbersome task. Therefore, the quantum entanglement properties involving an 
infinite chain structure is quite important, not only because the mathematical calculation is cumbersome but also because 
real materials are well represented by an infinite chain. Thus, in this paper we consider an entangled diamond chain with 
Ising and anisotropic Heisenberg (Ising-XXZ) coupling. Two interstitial particles are coupled through Heisenberg coupUng 
or simply two-qubit Heisenberg, which could be responsible for the emergence of entanglement. These two-qubit Heisenberg 
operators are interacted with two nodal Ising spins. An infinite diamond chain is organized by interstitial- interstitial and 
nodal-interstitial (dimer-monomer) site couplings. We are able to get the thermal average of the two-qubit operator, called 
the reduced two-qubit density operator. Since these density operators are spatially separated, we could obtain the concurrence 
(entanglement) directly in the thermodynamic limit. The thermal entanglement (concurrence) is constructed for different values 
of the anisotropic Heisenberg parameter, magnetic field and temperature. We also observed the threshold temperature via the 
parameter of anisotropy, Heisenberg and Ising interaction, external magnetic field, and temperature. 



I. INTRODUCTION 

Quantum entanglement is one of the most attractive 
types of correlations that can be shared only among 
quantum systems [1]. In recent years, many efforts have 
been devoted to characterize qualitatively and quantita- 
tively the entanglement properties of condensed matter 
systems, which are the natural candidate to apply for 
quantum communication and quantum information. In 
this sense, it is very important to study the entanglement 
of solid state systems such as spin chains [2] . The Heisen- 
berg chain is one of the simplest quantum systems, which 
exhibits entanglement; due to the Heisenberg interaction 
it is not localized in the spin system. 

In the last decade, several diamond chain structures 
have been discussed. It is interesting to consider the 
quantum antiferromagnetic Heisenberg model on a gener- 
alized diamond chain, which describes real materials such 
as Cu3(C03)2(OH)2, known as natural azurite. Honecker 
et al. [3] studied the dynamic and thermodynamic prop- 
erties for this model. In addition, Pereira et al. 0, [5| 
investigated the magnetization plateau of delocalized in- 
terstitial spins, as well as magnetocaloric effect in ki- 
netically frustrated diamond chains. More recently, Lis- 
nii [6*1 studied a distorted diamond Ising-Hubbard chain, 
and that the model also reveals the geometrical frustra- 
tion. Thermodynamics of the Ising-Heisenberg model on 
a diamond-like chain was widely discussed in the Refs 

The motivation for the study of the Ising and 
anisotropic Heisenberg (Ising-XXZ) model on a diamond 
chain, according to experiments of the natural mineral 
azurite, theoretical calculations of the Ising-XXZ model, 
as well as the experimental result of the dimer (interstitial 
sites) exchange parameter is based on a number of recent 
works. The 1/3 magnetization plateau and the double 
peaks both in the magnetic susceptibility and specific 
heat were observed in the experimental measurements 
prl , fT2| . It should be noticed that the dimer (intersti- 



tial sites) exchange is much stronger than those nodal 
sites. Various types of theoretical Heisenberg- model ap- 
proximate methods were proposed: renormalization of 
the density matrix renormalization group of the transfer 
matrix, density functional theory, high-temperature ex- 
pansion, variational mean-field-like treatment based on 
the Gibbs-Bogoliubov inequality, and Lanczos diagonal- 
ization on a diamond chain to explain the experimen- 
tal measurements (magnetization plateau and the double 
peaks) in the natural mineral azurite [13]. The localizable 
entanglement (LE) has been calculated for ground states 
of arbitrary Hamiltonians numerically by making use of 
the density renormalization group formalism and exhib- 
ited characteristic features at a quantum phase transi- 
tion. The LE is completely characterized by the maxi- 
mal connected correlation function for ground states of 
spin- 1/2 systems [14]. All of these theoretical studies 
are approximate. There is another possibility. Since 
dimer interaction is much higher than the rest, it can 
be represented as an exactly solvable Ising-Heisenberg 
model. In addition, experimental data on the magne- 
tization plateau coincide with the approximation Ising- 
Heisenberg model 0, 0, m, [Hi. 

Several studies have been done on the threshold tem- 
perature for the pairwise thermal entanglement in the 
Heisenberg model with a finite number of qubits. Ther- 
mal entanglement of the isotropic Heisenberg chain of 
spin has been studied in the absence Jl^and in the pres- 
ence of an external magnetic field [18|, |l9|. The pairwise 
thermal entanglement of nearest-neighbor qubits is inde- 
pendent of the sign of exchange constants and the sign of 
magnetic fields in the XX even-number qubit ring with 
a magnetic field. The thermal entanglement in Heisen- 
berg XX decreases with increasing temperature and the 
threshold temperature value Tth is independent of the 
external magnetic field. Although in some references the 
threshold temperature Tth is called the quantum critical 
temperature, here we call it the threshold temperature to 
avoid misunderstanding with the magnetic phase transi- 
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Figure 1: (Color online) Schematic representation of XXZ- 
Ising diamond chain. Red line represents the quantum bipar- 
tite coupling. 



tion temperature. 

It is quite relevant to study the thermal entanglement 
of the Ising-XXZ on diamond chain, as Ananikian et al. 

have discussed the thermal entanglement of the 
Ising-Heisenberg chain in the isotropic limit in the cluster 
approach. The paper is organized as follows: In Sec. II 
we present the Ising-XXZ model on diamond chain. We 
obtain the exact solution of the model via the transfer- 
matrix approach and its dimer (two-qubit) reduced den- 
sity operator in Sec. III. In Sec. IV, we discuss the 
thermal entanglement of the Heisenberg reduced density 
operator of the model, such as concurrence and thresh- 
old temperature. Finally, concluding remarks are given 
in Sec. V. 



II. ISING-XXZ DIAMOND HAMILTONIAN 

In this work, we study the Ising-XXZ model with 
mixed nodal Ising spins and interstitial anisotropic 
Heisenberg spins in the presence of an external mag- 
netic field on a diamond-like chain. The thermodynamic 
properties were previously discussed in Refs. [6-9]. A 
diamond-like chain is schematically illustrated in Fig 1. 
The Hamiltonian operator can be expressed as follows: 

N 



(1) 



where (^^,.,^6,^)^ = S^.Sl^ + Sl^Sl^ + ^Sl^Sl, cor- 
responds to the interstitial anisotropic Heisenberg spins 
coupling (J and A), while the nodal-interstitial (dimer- 
monomer) spins are representing by Ising-type exchanges 
(Ji). The Hamiltonian also includes a longitudinal exter- 
nal magnetic field acting on Heisenberg spins and a 
magnetic h acting on Ising spins. For convenience, we 
will consider the case ho = h. 

The quantum Heisenberg spin coupling can be ex- 



pressed using matrix notation. We have 
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We obtain the following eigenvalues by diagonalization 
of the two-qubit Heisenberg (dimer) exchange and by as- 
suming fixed values for //^ and /ii+i: 
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where their corresponding eigenstate in terms of standard 
basis {| tt), I II), I tl), I It)} are given respectively by 
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In quantum information the states \(p2) and \(ps) are 
known as the "magic basis" or Bell state, where two- 
qubit Heisenberg are maximally entangled. These states 
are responsible for the rise of entanglement at finite tem- 
perature. 

In order to study the two-qubit entanglement, we use 
the concurrence (2Q| , l2ll | , which is defined by 



C = K^i|^i)|, 



(9) 



where \ipi) means the spin- flip state. Therefore, the con- 
currence C will be for states \(p2) and l^^s), which corre- 
sponds to qubits that are maximally entangled. If the 
states \(pi) and {(p^) have no concurrence (C = 0), that 
means the qubits are unentangled. Later we give the 
definition of entanglement in more detail. 



A. Zero temperature phase diagram of entangled 
state 

In this section we study the phase diagram of the en- 
tangled states, similar to the magnetic phase previously 



2 





3 




2 


UFi 


1 


<^ ENT 







A 

-1 


UFM \. fa) 


-2 




-3 




-4 



ENO 










/ 


UFI 


UFM 


J 




(b) 





UFI ^ 








ENQ /'// 

/ / 
/ 

/ 


UFM 
(c) 



Figure 2: Phase diagrams of entangled state at zero temper- 
ature. The entangled state in frustrated phase denoted by 
ENT, whereas the entangled state in quantum ferrimagnetic 
phase we denote by ENQ. The unent angled state in ferri- 
magnetic phase (UFI) and unentangled ferromagnetic (UFM) 
states. 



discussed in Refs. [7','9*|, where it was observed three mag- 
netic states, a frustrated (FRU) state, a ferrimagnetic 
(FIM) state, and a ferromagnetic (FM) state. Although 
we only have two states (entangled and unentangled) , we 
discuss here the entangled region which is closely related 
to the magnetic states. It can be expressed the following 
states: 
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^3)1^ I and its eigenvalues is Eenq = -^-t~Jj' 
We also have the unentangled state UFI whose eigenvalue 
becomes Eufi = ~1 + T ~ ^^^^ phase boundary 
between UFM and UFI is given by h/J = 2, whereas 
between the ENQ and UFI states it is limited by A = 1, 
while for the ENQ and UFM states we have A = ^ - 3. 

Finally, in Fig. ^c) we illustrate the phase diagram 
Ji/J versus h/JdX zero temperature for A = 1 and 
h/ J = Hq/ where the eigenvalues are given by Eenq = 

3 h rp 1 Ji h rp liJi 3h 

-4 ~ 2J ' ^UFi - 4 - — - 2j and EuFM - 4 + " ~ 2J- 
The phase boundary between the ENQ and UFI states 
is limited by Ji/J= 1 , while the boundary between UFI 
and UFM is given by Ji = /i/2, and between the ENQ 



and UFM states the boundary follows the curve ^ 
J — 1. The doted line corresponds to the phase diagram 
for A = 0.5, whereas the dashed line represent the phase 
diagram for A = 1.5. 

However, entangled states will change with increasing 
temperature. This will be discussed in Sec. IV. 



III. THE PARTITION FUNCTION AND DEN- 
SITY OPERATOR 

With the aim of studying any thermal quantities, we 
first need to obtain a partition function on a diamond 
chain. As mentioned earlier (22l-[23|, this model can be 
solved exactly through a decoration transformation and 
transfer-matrix approach [26]. In order to summarize 
this approach we will define the following operator as a 
function of Ising spin particles /i and /i'. 



where | /i)^ stands for arbitrary values (/i = ±1/2) of the 
nodal spin in the i-th block. 

According to the Bell states given in eqs.® and 
([7j), at zero temperature the entangled (ENT) state is 
fully spanned in the frustrated state for zero magnetic 
field, while the ferrimagnetic and ferromagnetic state is 
spanned by unentangled states denoted by UFI and UFM 
respectively. 

In Fig. [2](a) we plot the phase diagram A versus Ji/J 
at zero temperature for h/J = ho /J = , where the 
ground-state energy of the entangled state is Eent = 
^ ^, while for the unentangled state (ferrimagnetic 



2 4 
state) it is Eufi 



J 



^, and for the unentangled 



state (ferromagnetic state) it is given by Eufm = ^ + 
The phase boundary between UFI and UFM is simply 
Ji/J = 0, and the phase boundary between the UFI and 
ENT states is given by (^ = ^ + ^), while the phase 
boundary between UFM and ENT (frustrated) state is 
given by (^ = -f - \)- 

In Fig. [2|(b) we plot the phase diagram A versus h/J 
at zero temperature for Ji/ J = 1 and h/ J = ho/ J^ where 
the eigenvalues for the unentangled state UFM is given by 
Eufm = 1 + t ~ §7' while for the entangled state (quan- 

N 

tum ferrimagnetic state) denoted by | ENQ) = W \ 



^(/i,/iO 



(13) 



where Habij^^j^') corresponds to the rth-block Hamilto- 
nian ([!]) (without summations), which depends on the 
neighboring Ising spins /i and /i^ Alternatively the oper- 
ator (p!3|) could be written in terms of two-qubit operator 
eigenvalues dU, which is 



g{p,,,') = Y,e-^'^^^'^'^^i){^i\. 



(14) 



Straightforwardly, we can obtain the Boltzmann factor 
by tracing out over the two-qubit operator. 



(15) 



where the Ising-XXZ diamond chain partition function 
can be written in terms of Boltzmann factors, 

Zn = ^^^;(/il,/i2) • . . ^(/iAr, /ii). (16) 
{m} 

Using the transfer-matrix notation, we can write the 
partition function of the diamond chain straightforwardly 
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by Zjv = tr (W^) , where the transfer-matrix is ex- 
pressed as 



W 



(17) 



where the transfer matrix elements are denoted by 
W-^-^ = = w{^^—^) and w = 

After performing the diagonalization of the transfer 
matrix ([l7j), the eigenvalues are 



W- 



±1 



(18) 



by tracing out over all Heisenberg spins and Ising spins 
except at the rth block in Heisenberg spins on a diamond 
chain. Using the transfer-matrix approach, we are able 
to write the reduced density operator by the following 
expression 



(22) 



Alternatively, using the transfer-matrix notation, we 
can write the reduced density operator as 



■ Aw\_. There- 



assuming that Q = y ('^++ — w 

fore, the partition function for finite chain under periodic 
boundary conditions is given by 



where we are assuming 
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Zn = 



A 



N 



(19) 



In the thermodynamic limit the partition function will 
be simplified, which results in Zn = A;?^. It is worth 
noticing that this result was widely used in Refs. [6-9]. 



A. Two-qubit operator 

In order to calculate the two-qubit Heisenberg opera- 
tor bonded by Ising particles /i and ja' ^ we assume that 
the Ising spin of the particle is fixed. Thus, the qubits 
operator elements in the natural basis becomes 



^1,1 

^2,2 92,3 

93,2 93,3 

94,4 



(20) 



by 



where the elements of the two-qubit operator are given 



^?4,4(m,m') =e-''^^(^'^'). 



(21) 



The thermal average for each two-qubit Heisenberg op- 
erator will be used to construct the reduced density op- 
erator. 
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The corresponding matrix U that diagonalizes the 
transfer matrix W can be given by 



U = 



A+ - 



and 



A_ - 



A_ —w- 
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(25) 
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Finally, the reduced density operator defined in eq. (|22]) 
must be expressed by 



tr ^U-^PijU 
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A^-i 
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(27) 



This result is valid for arbitrary number N of cells in a 
diamond chain, under periodic boundary conditions. 



C. Reduced density operator in thermodynamic 
limit 



Real systems are well represented in the thermody- 
namic limit {N 00), hence, the reduced density opera- 
tor elements after some algebraic manipulation becomes 



B. Reduced density operator and transfer matrix _ ^ / g^,j(2 > 2^^^^'J^~ 2 '~ 2-* _^ '^Qi,ji2 '~ 2^^^+- 

approach A+ \ ^ ^ 
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We will perform the reduced density operator bounded 2Q ( ' v^^) 

by Ising particles along a diamond chain. We can proceed ) 
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where we have assumed (A_/A+) ^ in thermody- 
namic hmit. 

Ah elements of reduced density operator immersed on 
a diamond chain are 

" 

p = ^' ^^'^ ^^'^ ^ . (29) 

^ P3,2 P3,3 ^ ^ 

_ P4,4_ 

It is worth noting that, this reduced density operator is 
the thermal average two-qubit Heisenberg operator, im- 
mersed in the diamond chain, and it can be verified that 
tr(/)) = 1. The cluster approach [15] becomes identical 
to transfer-matrix approach (the present approach) only 
when the magnetic field is zero. Later, in the next section 
(Fig. [8]) we show the difference between both approaches. 
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Figure 3: Density plot of concurrence C as function of T/J 
and A. Black color corresponds to the maximally entangled 
region, while by white color is the unentangled region, the 
gray region means the entanglement with degrees of different 
intensity In (a) we display C for absence of the magnetic field 
and in (b) we display C for h/ J — 1. 



IV. TWO-QUBIT 
MENT 



HEISENBERG ENTANGLE- 



Quantum entanglement is a special type of correlation, 
which only arises in quantum systems. Entanglement re- 
ffects nonlocal distributions between pairs of particles, 
even if they are removed and do not directly interact 
with each other. In order to measure the entanglement 
of anisotropic Heisenberg qubits in the Ising-Heisenberg 
model on a diamond chain, we study the concurrence (en- 
tanglement) of the two-qubits Heisenberg (dimer), which 
interacts with two nodal Ising spins using the definition 
proposed by Wooters et al. [20|, |2l|. 

The concurrence could be expressed in terms of a ma- 
trix 



(30) 



which is constructed as a function of the density operator 
p, given by Eq. ([29|) . with p*, which we represent by the 
complex conjugate of matrix p. 

Thereafter, the concurrence of two-qubits Heisenberg 
coupling (bipartite) could be obtained in terms of eigen- 
values of a positive Hermitian matrix R: 



C{p) = max{ V Al 



A4,0}, (31) 



with eigenvalues Ai ^ A2 ^ A3 ^ A4. Equation ([3T]) can 
be reduced to 



C{p) = 2max{p2,3 - ^/pi,ip4A, 0}. 



(32) 



A. Concurrence 



As presented above the entanglement can be studied in 
terms of concurrence defined by Eq. (|32]) . as a function 
of Hamiltonian parameters defined in Eq. ([T]) , as well as 
temperature and external magnetic field. 
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Figure 4: Density plot of concurrence C as function of h/ J and 
Ji/J . In (a) we display of A = 1.0 and (b) is for A = 2.0. 
The black (white) region corresponds to C = 1(0), and by 
gray regions we indicate a concurrence < C < 1. 



In Fig. [3] we illustrate the density plot of concur- 
rence C as a function of T/J and A for a fixed value 
of Ji/J=l. Black corresponds to the maximum entan- 
gled region (C = 1), while white is the unentangled region 
(C = 0). Gray means the different degrees of entangle- 
ment (0 < C < 1). After that we use this representation 
for the concurrence C, depending on the parameters of 
the Hamiltonian ([1]). The entangled region represented 
by different gray intensities depends also on the temper- 
ature. For high temperature the fuzzy region increases, 
while for the low-temperature phase between the entan- 
gled region and the unentangled region the boundary be- 
comes sharper. In Fig. it is shown that the model 
is maximally entangled only for A ^ 1 in the absence of 
the magnetic field, while the concurrence C is always less 
than 1 for A < 1 . The concurrence C becomes smaller 
with increasing temperature and the entanglement dis- 
appears at high temperatures. In Fig. [D^b) we display 
C for h/J = 1, where the concurrence behaves similar 
to Fig. [31(a). But there is the maximum entangled re- 
gion Fig. [31(b) for A ^ 1.0 and at temperatures less than 
T/J ^0.2. 

The density plot of concurrence C as a function of mag- 
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Figure 5: Density plot of concurrence C as a function of h/J 
and T/J: In (a) we display C for a A = 1 and in (b) we 
display C for A = 1.5. The black (white) region corresponds 
to C = 1(0), and by gray regions we indicate a concurrence 
< C < 1. 
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Figure 6: (Color online) Concurrence as a function of tem- 
perature for a fixed value of Ji/ J=l: In (a) we displays con- 
corrence for A = 2 and (b) displays concorrence for A = 0.9. 



netic field h/J and Ji/J is shown in FIGH] . To repre- 
sent the concurrence we use the same representation as 
in Fig. El for a fixed value of T/J = 0.3. In Fig. (Ka) 
we display C for a A = 1.0; thus, we can show that the 
concurrence is always less than C < 0.5. The largest en- 
tanglement occurs for small \Ji/J\ ^ 1 and small mag- 
netic field h/J < 2. Whereas in Fig. [J^b) we display 
C for A = 2.0, and the entanglement becomes stronger 
than for A = 1, but the concurrence is still limited to the 
regions \Ji/J\ ^ 1.5 and h/J<3. 

Another case is shown in Fig. [5] for the density plot 
of concurrence C as a function of h/J and T/J, at a 
fixed value of Ji/J = 1. Black (white) region indicates 
the entangled (unentangled) region, while the gray region 
corresponds to < C < 1. The concurrence is not very 
significant for A = 1, as shown in Fig. [5fa). For values of 
^ /i/J < 2 the entanglement becomes C < 0.25 at small 
temperatures, while at higher temperatures it disappears. 
For A = 1.5 the situation has changed; we observe the 
maximally entangled region for values of ^ h/J < 2.3 
and T/J< 0.2, as shown in Fig. 5 (b), and at higher 
temperatures and a strong magnetic field entanglement 
vanishes asymptotically. 

The concurrence as a function of temperature T/J for 
fixed values of Ji/J = 1 is shown in Fig. [6l for several 



Figure 7: (Color online) Concurrence as a function of 
anisotropic parameter A:(a) /i/J = and (b) h/J = 2.5 with 
different fixed values of temperature. 



values of the magnetic field. First of all, we illustrate 
the concurrence for A = 2 in Fig. [61(a). The concur- 
rence reaches its maximum value C = 1 in the absence 
of magnetic field and at low temperature, whereas at 
high temperature entanglement disappears, the threshold 
temperature occurs at T/J ^ 1.1, which is represented 
by a dash-dotted line. The behavior of the concurrence 
varies in a strong magnetic field, and the entanglement 
curve in most cases is lower than in the absence of a 
magnetic field. In the case of h/J = 2.4 (solid line), the 
maximum value of concurrence C = 1 is only achieved at 
very low temperature T/J, although the threshold tem- 
perature T/J is slightly larger than that in the absence 
of a magnetic field. In the case of the magnetic field 
be h/ J = 2.5 (dashed line), the concurrence achieves the 
highest value around C = 0.5 at zero temperature. While 
for h/ J = 2.6 (dotted line), the system becomes unentan- 
gled at zero temperature, but at low temperatures near 
T/ J ^ 0.3 it reaches its maximum value around C ^ 0.38. 
Second, in Fig. [61(b) the behavior of the concurrence for 
a fixed value of the anisotropic parameter A = 0.9, and 
for different magnetic field values are shown. Now, the 
behavior of concurrence decreases with different values 
of the magnetic field and has a peak near T/J ^ 0.3. 
The concurrence for all of these magnetic fields are un- 
entangled at zero temperature and there is a tempera- 
ture threshold near to Tth/J « 0.64 for ^ h/J < 2.4, 
while for h/J > 2.5 the threshold temperature jumps to 
Tth/J > 0.75. 

In Fig. [TJa) the concurrence as a function of the 
anisotropy parameter A in the absence of a magnetic 
field for various fixed temperatures is plotted. For A < 1 
the system is unentangled, and for A > 1 the system be- 
comes entangled. The maximum of the entangled region 
is achieved for larger anisotropic parameters, but when 
the temperature increases the entanglement decreases. 
The concurrence as a function of the anisotropy param- 
eter A, and for a fixed magnetic field h/J = 2.5 is plot- 
ted in Fig. [T^b), for the same set of temperatures as 
shown in Fig. [71(a). Here we can observe the influence of 
the magnetic field for the concurrence; at low tempera- 
ture (T/J = 0.01) the system is entangled for the values 
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Figure 8: (Color online) (a) Concurrence against exter- 
nal magnetic field, for different values of temperature with 
anisotropic parameter A = 1. (b) Concurrence as a function 
of temperature assuming fixed values of the anisotropic pa- 
rameter and a fixed magnetic field hJ — 0.2. By the solid 
lines we represent the cluster- approach results whereas 
by the dashed lines we represent our result. 



around of A « 2, whereas for higher temperature a small 
entanglement appears even for A < 1. 

Finally, in Fig. [HKa) we display the concurrence as a 
function magnetic field, assuming a fixed value of the 
anisotropy parameter A = 1. In this plot the differ- 
ence between the cluster approach ^| and the transfer- 
matrix approach used here is also shown. By the solid 
lines we represent the transfer-matrix approach, while 
with dashed lines we denote the cluster approach [15]. 
While, in Fig. [8Kb) the concurrence as a function of tem- 
perature (T/J) is illustrated, by the continuous line we 
observe the cluster approach, and by the dashed line we 
represent our current result. For the anisotropic param- 
eter A = 1.1 and for fixed magnetic field h/J = 0.2. 
Clearly, the cluster approach shows a double peak (solid 
line), whereas by dashed lines we represent our current 
approach discussed here. Therefore, we can conclude that 
the small peak is absorbed due to the Ising coupling of 
the diamond chain structure. However, the small peak is 
rather irrelevant since the concurrence maximum value 
is around C ^ 0.01. 

The difference between both approaches is consistent, 
since the entanglement on the diamond chain must be 
lower than for the cluster approach, where only the near- 
est Ising spins are considered, and we ignored the re- 
maining diamond blocks coupling contributions, so the 
cluster approach is equivalent to infinite uncoupled dia- 
mond blocks. 



B. Entanglement threshold temperature 

In Fig. [9](a) we display the phase diagram of the entan- 
gled region and the untangled region, as a function of the 
Ising coupling parameter Ji / J against threshold temper- 
ature for a fixed values of magnetic field h/J = 0.1, and 
for several anisotropy parameters. For each curve the 
left side (T/J < Tth/ J) of this region corresponds to the 
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Figure 9: (Color online) Phase diagram, of Ising coupling 
parameter Ji/J as a function threshold temperature Tth/J- 
(a) For several values of anisotropy parameter A. (b) For 
several values of external magnetic field h/J. 



entangled state, while for right side (T/J > Tth/J) the 
system becomes an unentangled region. The entangled 
region depends of the anisotropy parameter A; as long as 
the anisotropy parameter increases the entangled region 
increases, too. For a large negative value of Ji/J < 0, 
the unentangled region is limited around Tth/J ^ 0.3, 
whereas, for Tth/J ^ 0.3 we still have an entangled re- 
gion, although this entanglement leads to a weakly en- 
tangled region for large negative value of Ji/J, even for 
small a anisotropy parameter A. However, for large pos- 
itive value of Ji / J the system becomes a fully unentan- 
gled region, but for large values of anisotropy parameter 
A tries to enhance the entanglement, although it always 
becomes an unentangled region for Ji/J large enough. 
In Fig. [ni^b) we display the phase diagram of entangled 
region and untangled region as a function of the Ising cou- 
pling parameter Ji / J against threshold temperature for 
a fixed value of the anisotropy parameter A = 1 , and for 
several values of magnetic field. For null magnetic field 
the entangled region is limited by the dotted curve, (the 
upper and lower side of this curve are symmetric), but as 
soon as the magnetic field is turned on, this symmetry 
is broken even for a low magnetic field like h/J = 0.2 
(described by the dash-dotted line), the system will al- 
ways be entangled for Ji/J ^ —1.5 and Tth/J ^ 0.28. 
For a stronger magnetic field this limit increases until it 
achieves the asymptotic limit Tth/J ^ 0.910239. The 
largest value of threshold temperature Tth/J ^ 0.910239 
occurs when Ji/J = and for arbitrary values of mag- 
netic field. Of course, this limiting threshold temperature 
also depends on the magnitude of anisotropic parameters, 
as we can verify in Fig. [9Ka). 

Other important density plots, are displayed in Fig. 
[Tm for the concurrence C as a function of A against 
Ji/J, assuming a fixed temperature T/J = 0.3. This 
density plots follow same pattern (definition) as those 
illustrated in Fig. [3j the black (white) region corre- 
sponds to C = 1(0), while a gray region corresponds to 
< C < 1. In Fig. [TOlfa) the maximally entangled re- 
gion only appears for A > 1; we can also verify that 
the illustration is symmetric [(C(Ji/J) = C(— Ji/J))] 
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Figure 10: (Color online) Density plot concurrence C as func- 
tion of A versus Ji/J. In (a) we display C for a null magnetic 
field and in (b), we display C C for h/J — 1. 




Figure 11: (Color online) Density plot concurrence C as func- 
tion of T/J versus Ji/J. In (a) we display C for a /i/J = 
and in (b) we display C for h/J — 1. 



only when = 0. However, when the external mag- 
netic field is turned on [i.e., {h/J = 1)], the symmetry 
Ji/J ^ —J\/ J is broken as shown in Fig. [TOT b). When 
J\/ J > ^ the boundary between the entangled region and 
the unent angled region is sharper, while for Ji/J < the 
boundary between the entangled region and the unen- 
tangled region becomes less pronounced. The maximally 
entangled region becomes stronger for a large value of 
the anisotropy parameter A > 1 compared to the en- 
tanglement without magnetic field. This result is simi- 
lar to the previous result discussed by Zhou et al. [27], 
when they studied the anisotropy effect of a two-qubit 
Heisenberg model with external magnetic field, although 
they discussed only one-couple qubits. The red solid line 
is used to represent the boundary between the entan- 
gled region and the unentangled region. This curve is 
also known as the threshold-temperature curve; clearly 
for null magnetic field the threshold temperature rounds 
symmetrically the entangled region-this is not so evident 
when external magnetic field is switched on. Even for a 
rather weak external magnetic field, the entangled region 
spreads out for J\/ J < 0. The entanglement emerge for 
Ji/J < —1, [Fig. [TOlf b)]. despite the concurrence being 
tiny and at first glance the concurrence still looks like 
in Fig. [TOlfa). but the threshold temperature is highly 
different compared to that without magnetic field, such 
as there is no threshold temperature for Ji/J < — 1 and 
A > -0.3. 

The density plot of concurrence C as a function oiT / J 
versus Ji/J is depicted in Fig. [TTJ for a fixed value of 
A = 1.0. We show again that the maximally entangled 
region vanishes for higher temperature T/J ^ 0.3. In 
Fig. fTTTa) the density plot of concurrence C is displayed 
for /i/J = 0, and the maximally entangled region is lim- 
ited by I Ji/ J| < 1, following the same kind of density plot 
as in Fig. [TOl by the solid red line we represent the thresh- 
old temperature limiting the entangled region. Whereas 
for /i/ J = 1, as displayed in Fig. [TTT b). the maximally 
entangled region shrunk to —0.5 < J\/ J ^ 1, despite the 
entangled region spread out with weakly entangled region 
for Ji/ J < and T^^/ J < 0.5, consequently the symmetry 
in Ji/J ^ —Ji/ J is broken, this is better illustrated as 



a function of the threshold temperature curve described 
by the red solid line. 



V. CONCLUSIONS 

In general most of the quantum entanglement stud- 
ied of an infinite chain is a cumbersome issue, motivated 
by this fact we investigate the quantum entanglement of 
the Ising-XXZ diamond chain. Even when Ising coupling 
does not contribute directly in quantum entanglement, 
we discuss the Heisenberg two-qubit entanglement effect 
on Ising- Heisenberg diamond chain structure. First we 
obtain the Heisenberg dimer operators immersed in a di- 
amond chain with Ising coupling. Using this result we 
are able to obtain the average of the two-qubit operator. 
Thereafter, the concurrence is obtained straightforwardly 
in terms of the reduced density matrix operator elements. 
Using the concurrence, we study the entanglement of the 
Ising-XXZ diamond chain as a function of Hamiltonian 
parameters, such as temperature as well as external mag- 
netic field. 

Usually the entangled region vanishes when the tem- 
perature increases; in some cases [see for instance Figs. 
[TOlfb) and [TTT b)] the entanglement vanishes asymptot- 
ically when the external magnetic field is switched on. 
The entangled region is limited by the so-called thresh- 
old temperature, but for some other parameters [like in 
Figs. [TOlfb) and [TTTb)] the threshold temperature only 
occurs in the asymptotic limit. From our result we con- 
clude that there is no double peak in the concurrence as a 
function of temperature, such as that obtained by the use 
of the cluster approach[15] [see Fig. [8jb)]; the disappear- 
ance of the tiny peak can be understood as a dissipation 
due to the qubits are inmerse in diamond chain structure. 
However this small peak is irrelevant, since, the degree 
of entanglement is rather small: C = 0.01. 

It would be interesting also in the future to consider- 
ation, the case of tripartite entanglement for Heisenberg 
coupling of an Ising-Heisenberg chain instead of bipartite 
Heisenberg coupling. Somewhat similar to that consid- 
ered by Tsomokos et al. [28] where the tripartite coupling 
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was studied at the zero temperature. 
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